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Long-term polarimetric observations of OH127. 8+0.0 
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ABSTRACT 

OH 1612- and 1667-MHz masers from the well known object OH127. 8+0.0 were mon- 
itored in full-polarization mode over a period of 6.5 years and mapped with MERLIN 
at one epoch. The OH variability pattern of the star is typical of extremely long- 
period AGB stars. The distance determined from the 1612-MHz light curve and a new 
measurement of the angular radius is 3.87±0.28kpc. At both frequencies, the flux of 
polarized emission tightly follows the total flux variations while, the degrees of circu- 
lar and linear polarization are constant within measurement accuracy.. There is net 
polarization at both lines. The magnetic field strength estimated from a likely Zeeman 
pair is — 0.6mG at the distance of 5400 au from the star. At the near and far sides of 
the envelope, the polarization vectors are well aligned implying a regular structure of 
the magnetic field. The polarization characteristics of the OH maser emission suggest 
a radial magnetic field which is frozen in the stellar wind. 

Key words: masers - polarization - circumstellar matter - stars: AGB and post- AGB 



1 INTRODUCTION 

OH/IR stars have a spectral energy distribution which peaks 
between 10 and 25/im, an absorption feature at 9.7/im and 
a double-peaked OH maser profile at 1612 MHz. Almost all 
of them are long period variables with periods of 500-3000 
days and bolometric amplitudes of 1-2 magnitudes (Habing 
1996). As low and intermediate mass stars with enhanced 
mass loss up to a few 10~^Mq yr~^, which leads to a very 
optically thick and dusty circumstellar envelope, they could 
be considered as successors of Miras (van der Veen & Habing 
1988; Vassihadis & Wood 1993; Habing 1996). Untargeted 
surveys of OH 1612-MHz line suggested that OH/IR stars 
have progenitor masses from 1 to 8Mq (Baud et al. 1981). 
OH/IR objects with periods above 1000 days are thought 
to be probably medium mass (Js4.5Mq) evolved objects 
(Likkel 1989; Lepine, Ortiz & Epchtein 1995; Habing 1996). 
The study of mid-infrared spectral features in a large sam- 
ple of circumstellar OH masers provides indirect evidence for 
mass segregation between Miras and OH/IR objects (Chen 
et al. 2001). One of the ways to verify this hypothesis and to 
constrain the evolutionary status of OH/IR stars is to assess 
their parameters using reliable distances. 

OH127.8+0.0 (= OH127 = IRAS01304+6211 = 
V669 Cas) was discovered in a small-scale untargeted survey 
at 1612 MHz (Kerr & Bowers 1974). High angular resolution 
studies of this brightest OH/IR star in the anti-centre region 



E-mail: wolak@astro.uni.torun.pl 



revealed a thin OH 1612-MHz maser shell of ~l'.'5 radius 
expanding with a velocity of about 12kms~^ (Booth et al. 
1981; Norris, Diamond & Booth 1982; Bowers, Johnston & 
Spencer 1983). These observational characteristics are fully 
consistent with the predictions of the standard expanding 
shell model (Goldreich & Scoville 1976; Elitzur, Goldreich 
& Scoville 1976), although a large portion of the envelope 
shows departure from spherical symmetry. In this model, a 
thin spherical shell and double-peaked OH 1612-MHz pro- 
files are the natural consequence of radial amplification of 
OH maser emission in a uniform outflow at a constant veloc- 
ity (Reid et al. 1977). The star has received much attention 
over the years as an archetypal OH/IR object but its dis- 
tance is poorly determined spreading from 2.9 to 6.2 kpc 
(van Langevelde, van der Heiden & van Schooneveld 1990; 
Bowers & Johnston 1990). In this paper we estimate its 
distance comparing the linear diameter of the shell, deter- 
mined using the phase lag technique (Herman & Habing 
1985), with the angular diameter measured on interferomet- 
ric maps. 

In early observations, the OH maser emission from 
OH/IR stars appeared largely unpolarized (Cohen 1989). 
However, our recent study has revealed polarized features in 
about 60% (33/57) of the OH/IR stars (Wolak, Szymczak & 
Gerard 2012). EUiptically polarized emission with a degree 
of polarization up to 30% is usually detected. This indicates 
the presence of magnetic fields in the outer envelopes. High 
angular resolution observations have revealed a magnetic 
field of a few mG, whose structure was imaged in several 
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late-type stars (Chapman & Cohen 1986; Szymczak, Cohen 
& Richards 1998, 2001; Bains et al. 2003; Amiri, Vlemmings 
& van Langevelde 2010). OH127 has been mapped several 
times with the VLA (Bowers et al. 1983; Bowers & Johnston 
1990) and MERLIN (Booth et al. 1981; Norris ct al. 1982) 
but no full polarization data were obtained. We present po- 
larimetric maps of the target and analyze full polarization 
time series in order to explore the behaviour of the magnetic 
field and the degree of saturation of maser amplification. 



2 OBSERVATIONS 
2.1 Nangay 

Data were obtained at 189 epochs during the period 2002 
May 12 - 2008 December 28 (Table 1) with the Nangay 
radio telescope (NRT). A cooled 1.1— 1.8GHz receiver and 
a 8192-channel autocorrelation spectrometer were used for 
simultaneous observations of the two OH transitions at 
1612.231 and 1667.359 MHz. The receiver system temper- 
ature was about 35 K and the effective velocity resolution 
was 0.142 and 0.137 kms"^ at 1612 and 1667MHz, respec- 
tively. The half-power beam width was 3!5(RA) by 19'(Dec) 
and the conversion factor of the antenna temperature to the 
flux density for a point source was ~1.4K Jy~^. 

The two orthogonal linear polarizations and two oppo- 
site circular polarizations for each transitionwere measured 
in the 8-bank mode of the spectrometer. The Stokes pa- 
rameters /, Q and V were directly provided by the system, 
while the Stokes U was extracted by a horn rotation of 
45° . The four Stokes parameters were used to derive: the 
linearly polarized flux density, p = \/ Q"^ + U"^ , degree of 
circular polarization, the mc = degree of linear po- 

larization, mL = p/I and the polarization position angle, 
X = 0.5tan~^(f7/(5). The absolute and relative accuracies 
of gain measured with a noise diode were ~5% and ~1%, 
respectively. The error in the polarized intensity caused by 
the polarization leakage between the orthogonal feeds was 
about 2%. The baseline subtraction was done by frequency 
switching: this mode introduces an error of less than 0.6% 
in the polarization parameters, as compared to the position 
switching mode. W12 (an unpolarized source in absorption) 
and W30H (a strongly polarized maser source) were reg- 
ularly observed throughout the whole project in order to 
measure the instrumental polarization and consistency of 
the amplitude calibration. The absolute flux density scale 
was accurate to within 7—8%. A detailed description of the 
methods of observations, full polarization calibrations, anal- 
ysis of errors is given in Szymczak & Gerard (2004). A total 
integration time for the two horn positions of ~12min re- 
sulted in a 3f7rms sensitivity level in the Stokes I of ~0.4 Jy. 
All the velocities given in this paper are relative to the local 
standard of rest. 

Since the upgrade of the NRT in 2001 (van Driel, Pez- 
zani & Gerard 1996) interferences from local and/or satel- 
lite transmitters were reduced dramatically. Throughout the 
whole observing period, we noticed them at a few epochs 
only and the total loss of data was less than one per cent. 



2.2 MERLIN 

0H127 was observed in the 1612- and 1667-MHz OH tran- 
sitions on 2001 December 21 using the six telescopes of the 
MERLIN axray in full-polarization mode. The narrow-band 
(0.25 MHz) observations of the target switched at intervals 
of 3 minutes between the two maser lines and were inter- 
leaved with 1.5-min scans on 0141+579, the phase reference 
source, in wide-band mode (14 MHz) at the appropriate fre- 
quencies to obtain the optimum signal-to-noise ratio. The 
source 3C84 was observed both in narrow-band and wide- 
band modes in order to determine the phase offset correction 
that arises due to observations made in different frequency 
setups and to derive corrections for instrumental gain vari- 
ations across the bandpass. The flux density scale and the 
absolute polarization position angles were determined by ob- 
servations of 3C286. 

The data reduction weis carried out following standard 
procedures for MERLIN OH polarimetry (Diamond et al. 
2003) using the d-programs at Jodrell Bank and the As- 
tronomical Image Processing System (AIPS). The /, Q, U 
and V data cubes were obtained with a beam of 0'.'17x0'.'15. 
The spectral channel width was 0.5 kHz, which corresponds 
to 0.10 and 0.09 km s"^ at 1612 and 1667 MHz respectively. 
The l(Trms noise level in emission-free Stokes / single channel 
map was typically 17mJy beam~^. The error on the absolute 
flux scale is 5 per cent. The errors on the polarization angle 
were usually lower than ±2?0 for the strongly polarized com- 
ponents. The absolute maser position error was ±20 mas. 



3 RESULTS 

3.1 Total flux density curves 

The variations of the total integrated flux densities at 1612 
and 1667 MHz arc shown in Figure 1. In order to represent 
the light curve shape a sum of two periodic functions of the 
type discussed by David, Etoka & Le Squeren (1996) was 
fitted to the data. The time of maximum for the red-shifted 
integrated flux density. To, the period, P, the amplitude 
measured in the logarithmic scale. Am, and the asymme- 
try defined as the fraction of the period between minimum 
and maximum are given in Table 1. There are no signifi- 
cant differences in those parameters between the 1612 and 
1667-MHz lines with the exception of Am, which differs by 
a factor of 1.7. The ratios of the integrated flux at the two 
lines range from 5.4 at minimum light to 9.9 at maximum. 
This implies a partial suppression of the 1667-MHz emission 
as the brightness increases due to a competitive gain effect 
(Gray & Field 1996). 

The total flux density at the two frequencies shows 
smooth and regular variations with the exception of a 
bump on the rising branch of the light curve near MJD 
= JD-2450000 = 4450. A similar event occurred near JD 
= 2446800 as might be deduced from the 1612-MHz curve 
shown by van Langevelde et al. (1990) (their Fig. 2). The 
asymmetry deduced from our data is higher than that re- 
ported previously (Herman & Habing 1985; van Langevelde 
et al. 1990). The period inferred from our data is shorter 
by 434 days than reported in Herman & Habing (1985). 
It is also shorter by 38 days than the value given in van 
Langevelde et al. (1990) but is well within their uncertain- 
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Figure 1. Integrated flux density of tlie OH 1612- (upper panel) 
and 1667-MHz (lower panel) spectra of OH127 as a function of 
time. The uncertainties of the measurements are comparable with 
the symbol sizes. The blue line shows the model light curve. The 
arrow marks the date of the MERLIN polarimetric observation. 



Table 1. Measured and derived parameters for OH127. 8+0.0. 



Nartfay data 

Span of monitoring data (MJD) 2355 - 4735 
Number of observations 189 

Expansion velocity, 14(1612) (kms"^) 11.13 ±0.14 

Expansion velocity, 14(1667) (kms"!) 11.31 ±0.14 

Stellar velocity, V, (kms'^) -54.78 ± 0.14 

Time of maximum, To (1612) (MJD) 3164.5±2.5 

Time of maximum. To (1667) (MJD) 3086.1±20.2 
Period, P(1612) (days) 1599.8±4.1 
Period, P(1667) (days) 1596.1±4.8 

Amplitude, Am(1612) (mag) 1.416±0.011 

Amplitude, Am(1667) (mag) 0.829±0.013 

Asymmetry of the light curve (1612) 0.406±0.006 

Asymmetry of the light curve (1667) 0.343±0.013 

Phase lag, A*(1612) (days) 61.81±1.10 

SheU radius (au) 5358.3 ± 95.7 

MERLIN data 



ties. An ephemeris of the time of maximum based on Her- 
man & Habing and van Langevelde et al.'s data is inconsis- 
tent with that derived from our observations covering about 
one and half cycles. The large errors in their estimates can- 
not account for that inconsistency and we suggest signifi- 
cant changes of the hght curve of OH127 from one cycle to 
the next. We conclude that the OH maser curves of OH127 
strongly resemble the general characteristics of the optical 
light curve of long-period (<600 days) and extremely long- 
period (~1850 days) Miras (Lebzelter 2011; Groenewegen et 
al. 2009). 



RA offset of the shell centre'^ (mas) 
Dec offset of the shell centre" (mas) 
Outer shell radius, ro(1612) (mas) 
Inner shell radius, ri(1612) (mas) 
Expansion velocity, Vc(1612) (kms 
Outer shell radius. To (1667) (mas) 
Inner shell radius, ri(1667) (mas) 
Expansion velocity, Vc(1667) (kms 
Stellar velocity, Vt (km s ~ ^ ) 
Distance (kpc) 



0.0 
-191.0 
1380 ± 40 
1150 ±40 
11.7±0.2 
1340 ± 70 
1120 ± 70 
11.4 ±0.2 
-54.95 ±0.20 
3.87 ±0.28 



"Relative to the phase centre of 
RA(J2000) = 0l'^33™51?232, Dcc(J2000) 



62''26'53'.'236 



3.2 Phase lag 

The time shifts between the 1612-MHz flux curves of individ- 
ual channels, i.e. phase lags are determined using the follow- 
ing method. The light curves of each channel were normal- 
ized and smoothed with a three-point running mean and lin- 
early interpolated to obtain uniformly one-day spaced data. 
The channel with the strongest emission at — 65.9kms~^ 
was chosen as a reference, i.e. as the zero point of the phase 
lag. All the other channels were scaled in amplitude and 
shifted in time to minimize the least squares difference be- 
tween their light curves and that of the reference channel. 
The phase lags of all channels are plotted in Figure 2. The 
linear relation between the phase lag and velocity is fully 
consistent with the model of a uniformly expanding shell. 
The phase lag between the strongest blue- and red-shifted 
peaks, A$, is 61.81±1.10 days (Table 1). Although the light 
curve is not strictly the same for each stellar cycle, the in- 
terpolation and smoothing applied to the data do not affect 
the phase lag because our observations were highly sampled 
with a typical time resolution of two weeks (~0.009P). 

Our value of A$ differs significantly from that of 
75.6±3.6 days and 41.4±8.1 days reported by Herman & 
Habing (1985) and by van Langevelde et al. (1990), respec- 
tively. Such a drastic difference between these previous esti- 
mates of the phase lag is likely to be due to their different 
methods as, the two teams used a similar data sets. Our 
method is similar to that of van Langevelde et al., but the 



data have a much higher temporal sampling. We confirm an 
asymmetry of the OH shell that may cause an error in the 
phase lag measurements of up to 15% (Bowers & Johnston 
1990). 

3.3 Polarization spectra and light curve 

Figures 3 and 4 show the average polarization spectra near 
the maximum and minimum and the time series of polar- 
ization parameters for the strongest features at 1612 and 
1667 MHz, respectively. 

The expansion velocity defined as the half-maximum 
width in the OH profile peaks at 1612-MHz, 14(1612) does 
not differ from that at 1667-MHz, Vo(1667) (Table 1). The 
stellar velocity determined as the mean velocity of the outer 
1612-MHz peaks, 14 (Table 1) and the OH expansion veloc- 
ities are quite consistent with previous estimates (Bowers & 
Johnston 1990). The outflow velocity of 13.0 km s~^ deter- 
mined from the observations of CO transitions (De Beck et 
al. 2010) suggests that the OH masers form in regions which 
still did not reach the terminal velocity. 

The 1612-MHz profile shows strong left hand circularly 
(LHC) polarized emission over almost the whole velocity 
range. The degree of circular polarization is usually lower 
than 20%. It generally decreases with an increase of the total 
flux density presumably due to depolarization effect. The 
degree of linear polarization is generally lower than 6% and 
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Figure 2. Phase lag between the light curves of spectral channels 
in the 1612-MHz spectrum of OH127 plotted against the velocity. 
The dashed black line is the best least squares fit to the whole 
data set. The red lines cross the best fit line at the velocities of 
the strongest blue- and red-shifted peaks. 

shows a slight difference between the blue- and red-shifted 
parts of the spectrum. The polarization position angle is 
very similar in all spectral channels with a mean value of 
-14?7±1?5. 

The 1667-MHz profile also shows strong LHC emission, 
but a relatively narrow polarized feature is seen only in 
the red-shifted part of the spectrum and rac is lower than 
12%. A linearly polarized feature is also detected in the red- 
shifted edge of the emission. The maximum value of ttil is 
6-7%, while the x angle changes rapidly across the profile of 
l.lkms^^ width, the mean value being 10?2±7?6. We con- 
clude that the quantitative polarization properties of the 
1612- and 1667-MHz maser profiles of OH127 are typical for 
OH/IR objects in the sample recently studied (Wolak et al. 
2012). 

The timeseries of all polarization parameters of the 
brightest 1612- and 1667-MHz features at the blue- and red- 
shifted parts of the spectrum (Figs. 3 and 4) exhibit regular 
variations of the Stokes V and p fluxes tightly correlated 
with the total flux variations. The mean and median val- 
ues of the degrees of circular and linear polarization and the 
polarization angle of those features are summarized in Ta- 
ble 2. The degrees of circular and linear polarization were 
constant within the uncertainty of instrumental polariza- 
tion of about 2-3%. The polarization angle also does not 
show any systematic variations higher than 8-13°. Outliers 
in the plots are of instrumental origin, probably due to un- 
balanced gains of both types of polarization that we could 
not correct using the calibration data available. We found 
that the difference between the average polarization position 
angles for the four brightest blue- and red-shifted channels 
of the 1612-MHz proflle was 1?0±13?5. This suggests that 
the 1612-MHz maser regions are permeated with a regular 
magnetic field that is stable over a ~6.5 year period. The 
linearly polarized emission at 1667 MHz was detected only 
in the red-shifted part of the profile. In the velocity range 
of -43.9 to -43.6 km s~^ the mean difference of x between 
1667- and 1612-MHz spectra is only 6?5±15?3 but it in- 
creases to 46?9±18?4 for the range -44.6 to -44.0kms~\ 
This indicates nearly the same projected direction of mag- 



Table 2 . Mean and median values of the polarization parameters 
for the brightest OH features derived from the multi-epoch data. 



Line 
(MHz) 


Velocity 
(kms-i) 


Parameter 


Mean 




Median 


1612 


-65.92 


mc(%) 


-5.34 


1.81 


-5.32 






mL(%) 


5.13 


1.90 


5.00 






x(°) 


-13.53 


6.97 


-25.02 




-43.65 


mc(%) 


-1.46 


1.68 


-1.41 






mL(%) 


2.53 


1.12 


2.45 






x(°) 


-23.32 


13.14 


-25.03 


1667 


-44.02 


mc(%) 


-7.20 


2.18 


-7.07 






mL(%) 


6.05 


1.56 


5.81 






x(°) 


12.95 


12.86 


11.82 



" Standard deviation 



netic fields in the compact regions of the far sides of the 
1612-MHz and 1667-MHz shells, whereas in the annular re- 
gions, these transitions likely probe different clouds. 

3.4 Shell parameters and distance 

The maps of OH emission at 1612 and 1667 MHz are well 
known (Booth et al. 1981; Norris et al. 1982; Bowers et al. 
1983; Diamond et al. 1985; Bowers & Johnston 1990) and 
we added the Stokes I spot maps (Fig. 5) only for compar- 
ison purposes. The new observation generally confirms the 
paucity of 1667-MHz emission in the eastern side of the well 
pronounced 1612-MHz shell. The extent of 1667-MHz emis- 
sion in the western side is the same as that of 1612-MHz 
shell, which indicates that the two masers operate at similar 
distances from the central star. 

The locations and velocities of the Stokes / maser com- 
ponents at both frequencies obtained from the MERLIN 
observations are used to refine the shell parameters: stel- 
lar position, expansion and stellar velocities, outer, ro and 
inner, ri shell radius. The standard thin-shell model (Reid 
et al. 1977) was fitted to the Stokes I spots (Fig. 6) using 
a weighted least-squares algorithm. Parameters for the best 
fitting model are listed in Table 1. These are consistent with 
those reported by Bowers & Johnston (1990). 

Assuming a spherical thin shell model with isotropic 
emission, the measurement of the 1612-MHz emission extent 
to the 3(7 level integrated over the velocity range from —56.0 
to — 54.2kms~^ reveals an OH shell radius of 1343±30mas, 
which agrees with the estimate of the standard model (Table 
1). The assumption of a spherical thin shell with isotropic 
emission allows us to consider this radius as the angular 
size of the shell. For non-isotropic emission, the angular 
size would be overestimated by ^10%. An increase of shell 
thickness and maser saturation leads to higher biases (van 
Langevelde & Spaans 1993). Combining the shell radius with 
the phase lag provided by the NRT observations, we obtain 
a distance of 3.87±0.28 kpc for OH127. 

3.5 Spatial structure of the polarized emission 

The polarization data obtained with MERLIN are analyzed 
using the standard methods (e.g Szymczak et al. 1998; Bains 
et al. 2003). Tables 3 and 4 list the parameters of the polar- 
ized components at both frequencies. Only the components 




Figure 3. Polarization spectra of OH127 at 1612 MHz. (a) The dashed and sohd hnes are the average spectra for MJD ranges 3072—3203 
(27 observations) and 4013—4168 (9 observations) near the maximum and minimum of the hght curve, respectively. The Stokes I, V, 
degree of circular polarization, mc, linearly polarized flux density, p, degree of linear polarization, rriL, and polarization position angle, 
X, are shown from top to bottom, (b) Time-series of the polarimetric parameters for the extreme blue- (—65.92 km s~^) and red-shifted 
(-43.65 kms-l) emission peaks. For the Stokes /, V and p the left- and right-hand ordinates give the flux density scales for the extreme 
blue- and red-shifted peaks, respectively. Error bars are shown only for the red-shifted emission peak for sake of clarity but those for 
the blue-shifted emission peak are similar. For the Stokes I values the error bars are generally smaller than or of comparable size to the 
symbols. 



with V and p brightness higher than ~ 100 mjy beam" ^ are 
listed. These make up about 10 and 7% of all Stokes / com- 
ponents at 1612 and 1667 MHz respectively. 

The spatial distribution of the 1612-MHz polarized 
components is shown in Fig. 7. The origin of the maps is 
the position of the central star as derived from the stan- 
dard model (Sect. 3.4). All of the 28 detected compo- 
nents have velocities larger than 0.65—0.7014 with respect 
to the stellar velocity and their typical Stokes V flux is 
0.1-0.2 Jybeam"\ There is a strong excess oi —V com- 
ponents (20/28). The blue-shifted components are much 
closer to the projected position of the central star than the 
red-shifted components. The degree of circular polarization 
ranges from 5 to 38% which is a factor of 2 higher than 
that observed in the NRT spectra implying depolarization 
in the NRT beam. Diflterence in the spectral resolution by 
a factor of 1.5 between both instruments also plays a role. 
There are only a few components with elliptical polariza- 
tion. The linearly polarized emission, generally weaker than 
0.5 Jy beam" ^ is detected from the blue-shifted part of the 



spectrum. The components with mL ^9% appear in a cluster 
of size ~400 mas located '-^200 mas west of the stellar posi- 
tion. The polarization angles range from —37 to 40°with a 
mean flux-weighted value of — 4?6±12?4 that is consistent 
with the NRT resuhs (Sect. 3.3). 

The polarized 1667-MHz components are seen at the 
red-shifted velocities. The circularly polarized components 
of flux lower than 0.3Jybeam~^ and \mc\ ^25% are clus- 
tered near the stellar position (Fig. 8). All of them appear 
as LHC components which is consistent with the NRT re- 
sults. The linearly polarized components show a similar loca- 
tion as the circularly polarized emission. The mean value of 
mL=23.6±5.3% is a factor 5 higher than that for the 1612- 
MHz emission. The mean flux-weighted polarization angle 
is -9?5±11?5. This agrees weU with that at 1612MHz and 
implies the same projected direction of magnetic fields in 
the regions of both OH masers. 

While the degrees of circular polarization at 1612 and 
1667 MHz are similar, the degrees of linear polarization at 
1667 MHz are much higher than at 1612 MHz. This can im- 
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Figure 4. Same as in Fig. 3 but for 1667-MHz line and with two exceptions: the average spectra are for MJD ranges 3081—3210 (27 
observations) and 4042—4224 (9 observations) near the maximum and minimum of the light curve, respectively and that the time-series 
of all polarimetric parameters are shown for the extreme red-shifted (—44.02 km s~^) emission peak, while the Stokes / flux density only 
is shown for the extreme blue-shifted (— 65.49 kms~^) emission peak. 



ply that the linear polarization at 1612 MHz is suppressed 
or depolarized in the regions of strong amplification of this 
line (Elitzur 1996). 



3. 5. 1 Magnetic field strength 

Two groups of spatially coincident (<0'.'l) RHC and LHC 
features at RA(J2000) = 0l''33™51';2247, Dec(J2000) = 
62°26'53'.'311 appear in 8 contiguous channels. These are 
marked by z in Table 3. The features are clearly asymmet- 
ric but can be easily fitted with two slightly shifted Gaus- 
sian components (Fig. 9). The strongest component near 
— 65.18 kms"'^ shows a 0.142±0.045 kms"'^ difference in the 
peak velocities between the RHC and LHC emissions. If one 
can interpret this difference as a Zeeman effect, then the 
magnetic field strength is — 0.60±0.19mG assuming an av- 
erage splitting coefficient of the a components (Davies 1974) . 
Here, the minus sign indicates that the direction of the field 
is towards the Earth. No linearly polarized emission is de- 
tected at this position. 



3.6 Emission from the near and far caps 

The 1612-MHz emission in the velocity ranges from —66.8 to 
— 66.0kms~^ (focap) and from —43.6 to — 42.8kms~^ (»"cap) is 
unresolved or barely resolved with the MERLIN beam. The 
lower limits to the brightness temperature (Tb) for &cap and 
^cap regions are 0.7—3.1x10^ and 3.3— 6.3 x 10^ K, respec- 
tively. Tb of the emission from fccap is a factor of 2— 7 higher 
than that of reap, while the angular size is by 18—23% lower. 
Both extreme regions coincide within 70 mas supporting the 
standard model of spherically symmetric shell. The emission 
from the near cap may contain a component from the am- 
plified stellar image (Norris et al. 1984; van Langevelde & 
Spaans 1993; van Langevelde et al. 2000). We failed to ver- 
ify this hypothesis due to the insufficient angular resolution 
of our data and the poorly known optical position of the 
underlying star. 

4 DISCUSSION 

4.1 Revised properties of OH127 

The present observations allow us to refine the basic stellar 
parameters of OH127. Whitelock, Feast & Catchpole (1991) 
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Figure 5. Stokes / maps of the 1612- and 1667-MHz maser spots 
of OH127. The symbol sizes are proportional to the logarithm of 
the brightness. The velocity scale is coded in colour. 
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Figure 6. Angular distance of the maser emission spots from the 
estimated stellar position of OH127 versus their radial velocity. 
The symbol sizes are proportional to the logarithm of the bright- 
ness. The least-squares fits of the thin shell model to the data are 
used to estimate the shell parameters as listed in Table 1. The 
red-dashed and blue-dotted curves show the outer radii fitted to 
the maser distribution at 1612 and 1667 MHz, respectively. 



Figure 7. Distribution of the 1612-MHz maser polarized compo- 
nents of OH127. Velocities are coded by colours as shown in the 
bottom right wedge. Upper panel: linearly polarized components. 
The size of squares is proportional to the logarithm of the p flux 
density and the length of polarization vectors is proportional to 
the degree of linear polarization, ttil. Lower panel: circularly po- 
larized components. The size of the circles is proportional to the 
logarithm of the Stokes V flux density. Empty and solid circles 
indicate negative and positive values of the Stokes V, respec- 
tively.The size of sectors is proportional to the degree of circular 
polarization. 



derived a period- luminosity, P — relation for OH/IR 
stars with pulsation periods of 1000—2000 days. Using their 
eq. 2, one would predict that the luminosity of OH127 is 
2.68li[ Qg X 10*1/0. The uncertainty is rather large but this 
value is close to 1.5 x 10'*Lq, which is the value adopted by 
De Beck et al. (2010) scaled to our distance estimate. Kem- 
per et al. (2002) obtained the best fit of the spectral energy 
distribution (SED) of OH127 to the observations at the dis- 
tance of 1.8 kpc for the assumed luminosity 0.63 x IO^Lq 
which translates to 2.9 x IO^Lq for our distance estimate 
(Tab. l).Thus, the stellar luminosity of OH127 is below the 
AGB-limit of 5.5 x lO" Lq (Wood et al. 1983) but it is very 
luminous compared to a typical AGE star. 

OH127 emits virtually all its radiation at infrared 
wavelengths and the mass- loss rate, M, obtained from 
modelling of the spectral energy distribution (SED) is 
1.5-2.8xlO"''MQyr"' (Justtanont & Tielens 1992; Kem- 
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Figure 8. Same as Fig. 7 but for the 1667-MHz emission. 
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Figure 9. MERLIN 1612-MHz spectra of the left (red solid line) 
and right (blue solid line) circular polarizations of maser fea- 
tures of close spatial coincidence. The dashed lines show the 
Gaussian curves fitted to the data. The strongest component 
near — 65.18kms~^ exhibits a velocity shift between the left (red 
dashed line) and right (blue dashed line) circular polarizations, 
which is a signature of Zeeman splitting. 



per et al. 2002) when scaled to our distance estimate and for 
an assumed dust to gas ratio of 0.01. The observations and 
modelling of multiple low- and high— J CO lines (up to J = 
7-6) by De Beck et al. (2010) give M = 3.4xlO"^M0 yr"^ 
for the distance of 3.87 kpc. The mass-loss rate calculated 
from an empirical formula based on the photo-dissociation 
model (Netzer & Knapp 1987) is 7.1x10"^Mq yr"^ This 
value is only a factor of two higher than that obtained with 
CO observations (De Beck et al. 2010) and is typical for 
OH/IR stars. This indicates that the M estimates from the 
SED modelling are an order of magnitude higher than those 
from CO and OH data, likely due to the poorly constrained 
dust to gas ratio (Justtanont & Tielens 1992). 



4.2 Magnetic field 

The strength of the magnetic field tentatively detected in 
OH127 is comparable in magnitude with that estimated for 
other stellar OH maser sources (Szymczak et al. 1998; Bains 
et al. 2003;Amiri et al. (2010); Etoka & Diamond 2010). 
However, it is an order of magnitude higher than that mea- 
sured from single dish high resolution spectra of several 
OH/IR objects (Zell & Fix 1991). This discrepancy is likely 
due to spatial blending of individual maser clumps in the 
envelopes and is confirmed for OH127: our NRT 1612-MHz 
spectra in the velocity range of —66.6 to — 64.5kms~^ con- 
tain features with mc ^5% while the MERLIN maps show 
highly (up to 38%) polarized components (Tab. 3). 

The presence of magnetic field in the outer envelope of 
OH127 is implied by the net polarization observed in both 
OH transitions. Specifically, the source shows a strong excess 
of LHC emission well seen in the NRT spectra and MER- 
LIN maps (Figs. 3 — 8). Enhancing one sense of circular 
polarization and inhibiting the other can be explained by 
the overlap of different Zeeman components due to a veloc- 
ity gradient within the medium (Deguchi & Watson 1986) 
or matching the gradients of the magnetic field and velocity 
(Cook 1966). The latter mechanism may be less probable 
in a thin circumstellar shell. A simple model of dipole mag- 
netic field with an adapted Cook's mechanism was successful 
in explaining the spatial segregation of the 1612-MHz spots 
dominated by the two senses of circular polarization in the 
supergiant star VXSgr (Zell & Fix 1996). It is quite pos- 
sible that Cook's mechanism could also partly explain the 
distribution of the 1612-MHz maser polarized components 
of OH127 (see Fig. 7) by adopting a scenario inspired from 
Zell & Fix's Fig. 4 but with dipole field lines stretched in 
the radial direction and a different inclination of the dipole 
axis to the line of sight. At a distance of 5350 an from the 
star a radial magnetic field or a field with open lines can oc- 
cur. In the first case the outfiow velocity gradient can com- 
bine with the Zeeman pair to produce a coherence length 
of maser amplification longer for LHC component than for 
RHC component at the near and far sides of the envelope, 
i.e. the blue- and red-shifted parts of the spectrum are dom- 
inated by LHC emission. If the field with open lines or of 
interstellar origin is oriented towards the observer, as shown 
by our observations, then the masers from the near and far 
sides of the envelope will be dominated by LHC and RHC 
emission, respectively. This scenario is not fully consistent 
with our data where the red wing shows positive and nega- 
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tive values of mc at 1612 MHz and only negative values at 
1667 MHz. 

OH127 at 1612 MHz shows very small differences in the 
polarization angles at the near and far sides of the shell, 
suggesting a regular configuration of the magnetic field. Such 
a phenomenon was observed in other OH/IR objects and 
two explanations are proposed in Wolak et al. (2012): (1) 
the field is of stellar origin and the same orientation of the 
field is carried away by the wind fiowing towards the near 
and back sides of the shell, (2) the field is of interstellar 
origin and can be amplified and/or distorted by the stellar 
wind. 

In the first case the field of 0.6 mG strength of solar 
type (B ~ r^^) gives a surface field strength of ~lkG for 
the assumed star radius of 859Rq (De Beck et al. 2010), 
which is ^ 10^ times stronger than the average magnetic 
field (~few G) for an AGB star (Kemball & Diamond 1991). 
In the second case, if the magnetic field pressure equals the 
gas kinetic pressure, then B^/Stt ~ pV^ , where p is the gas 
density in gcni^'^ Since p — M /A-Kr^Ve for a free radial 
expansion then B(G) ~ (2MV;)°-^/r. With M in units of 
10~®MQyr~^, Ve in units of lOkms"^ and the distance to 
the star, r in au we obtain 

S(G) = 2.4 X (2MV;)° V^- (1) 

For M = 3.4 (Dc Beck et al. 2010), Vc = 1.12, r=5350au 
(Table 1) this equation yields a B field strength of 0.9 mG, 
which is close to our value of 0.6 mG. This means that in 
our measured OH maser region there is nearly equilibrium 
between the magnetic and kinematic pressure. Eq. 1 may 
not be true close to the star where a solar (~ r~^) or dipole 
(~ r~^) type field may exist. On the other hand, the mag- 
netic field measurements from maser observations indicate 
that the field strength follows the relation B ~ n°'® over 
an enormous range of the number density, n, (Vlemmings 
2012). This empirical relation predicts a number density of 
molecular hydrogen of 1 x 10^ cm~^ for a 0.6 mG field. For 
the kinematic model (Sun & Kwok 1987) such a density is 
sufficient to sustain the 1612-MHz maser operation when 
M > 5 X 10"^ and > 1 x 1O~*M0 yr"^ for relative abun- 
dances of OH molecule higher than 2 x lO"'' and 5 x 10~*, 
respectively. These mass-loss rate values are consistent with 
those recently obtained with the CO observations and mod- 
elling (Dc Beck et al. 2010 and Sect. 4.1). 

We note that a magnetic field of 0.6 mG is two orders of 
magnitude larger than the average interstellar magnetic field 
(Troland & Heilcs 1986). This suggests that the magnetic 
field in the OH regions could be compressed or amplified by 
a factor of ~100 (Sokcr 2002). A model of field compression 
by Gerard (1985) predicts an amplification factor of up to 10. 
Another possibility is local enhancements of the magnetic 
field related to the stellar spots or convective cells which 
might produce the magnetic clouds similarly as in the solar 
wind (Soker 2006). 



4.3 Saturation level 

Our observations show that the lower limit to the bright- 
ness temperature of the strongest 1612-MHz emission is up 
to 6xl0^K (Sect. 3.6). VLBI observations at a phase of 
0.67 detected only the extreme blue-shifted emission near 



— 66.0kms^^, which was unresolved with ~100mas resolu- 
tion, while the extreme red-shifted emission was resolved 
out (Norris et al. 1984). Tb for the unresolved emission was 
estimated at 5 x 10® K. Our MERLIN data were obtained 
at phase 0.44, so it is likely that the flux density level was 
similar to that observed with the VLBI. Using the VLBI es- 
timates of size and brightness temperature for the extreme 
blue-shifted emission and the distance and shell radius given 
in Table 1 we estimated that the maser beaming solid angle, 
f7, is 0.016 sr. This is close to the canonical value (Goldre- 
ich & Keeley 1972). From our comparison of the angular 
sizes of the two extreme maser regions (Sect. 3.6) we get 

0. 022—0.024 sr for the beaming solid angle of the red-shifted 
emission. In the standard maser theory (Goldreich, Keeley 
& Kwan 1973) the 1612 MHz maser stimulated emission rate 
can be approximated (Amiri et al. 2010)) by 

«-'-'^'»-"(i^)(iA;)^" P' 

This implies R = 0.011 s"^ and 0.008-0.002 s"^ for the 
blue- and red-shifted emission regions, respectively. We no- 
tice that the first value is consistent with R predicted for 
interstellar OH masers (Goldreich et al. 1973). The radia- 
tive and collisional decay rate for the 1612 MHz is not obser- 
vationally constrained but its upper limit can be estimated 
(Elitzur 1992) with the following equation 

With the above estimated values of Tb and fl for the blue- 
shifted emission F < 1.5 x 10^'^ s^^ which implies R/T >7, 

1. e. the maser in OH127 is weakly saturated. This result ap- 
pears to be consistent with previous observations showing a 
fixed relation between the OH and infrared emissions (Har- 
vey 1974). The lack of relatively large random variations in 
the OH maser fluxes of OH127 also suggests that the maser 
emission is saturated. 



5 SUMMARY AND CONCLUSIONS 

Full polarization monitoring of 1612- and 1667-MHz OH 
maser lines over a period of 6.5 years and one epoch po- 
larimetric map of OH127 led to the following new results: 

(1) The OH maser fluxes show large (0.8— 1.4 mag) am- 
plitude variations with period of ~ 1600 days. The ratios of 
the integrated flux at 1612 and 1667 MHz lines of 5.4 and 
9.9 at the minimum and maximum, respectively indicate a 
partial suppression of the 1667 MHz emission. 

(2) The comparison of the light-travel diameter of the 
1612-MHz maser shell obtained using the NRT with the an- 
gular diameter determined from MERLIN maps yields a dis- 
tance of 3.87±0.28 kpc. 

(3) The fluxes of polarized emission follow the variations 
of the total flux, whereas the degrees of circular and linear 
polarization are constant within measurement accuracy. 

(4) The tentative detection of a Zeeman pair at 
1612 MHz implies a magnetic field strength of — 0.6mG at a 
distance of ~5400au. The strong alignment of polarization 
position angles at both the near and far sides of the shell and 
the presence of a net polarization suggest a regular magnetic 
field. 

(5) The stellar luminosity of OH127 is slightly below 
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the AGB limit and the mass loss rate is typical for an AGB 
star. 

(6) The measured sizes of the compact emissions from 
the near and far sides of the shell and the lower limits to the 
brightness temperature imply weakly saturated emission. 
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Table 3. Highly circularly polarized components in OH127. 8+0.0. The columns are as follows: (1) RA and DEC relative to the shell 
centre (Table 1); (2) the peak Stokes / flux density, together with the rms (Ic); (3) the peaJj Stokes V flux density, together with the 
rms (If); (4) the degree of circular polarization, together with the rms (1<t); (5) the peak p flux density, together with the rms (Icr); (6) 
the degree of linear polarization, together with the rms (Ic): (7) the LSR velocity of the component; (8) Zeeman pair. 
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Table 4. Highly linearly polarized components in OH127.8+0.0. The columns are as follows: (1) RA and DEC relative to the shell centre 
(Table 1); (2) the peak Stokes I flux density together with the rms (Ic); (3) the p peak flux density together with the rms (1<t); (4) the 
degree of linear polarization together with the rms (Icr); (5) the polarization position angle together with the rms (lu); (6) the peaJc 
Stokes V flux density together with the rms (Icr); (7) the degree of circular polarization together with the rms (If); (8) the LSR velocity 
of the component. 
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